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Abstract 

Exclusive semileptonic B decays into excited charmed mesons are investigated 
at order Aqcd /mg in the heavy quark effective theory. Differential decay 
rates for each helicity state of the four lightest excited D mesons {Di, D2, 
Dq, and Dl) are examined. At zero recoil, Aqqb/^q corrections to the matrix 
elements of the weak currents can be written in terms of the leading Isgur- 
Wise functions for the corresponding transition and meson mass splittings. A 
model independent prediction is found for the slope parameter of the decay 
rate into helicity zero Di at zero recoil. The differential decay rates are 
predicted, including Aqqb/^q corrections with some model dependence away 
from zero recoil and including order corrections. Ratios of various exclusive 
branching ratios are computed. Matrix elements of the weak currents between 
B mesons and other excited charmed mesons are discussed at zero recoil 
to order Aqcd/"Iq- These amplitudes vanish at leading order, and can be 
written at order Aqcd/w-q in terms of local matrix elements. Applications 
to B decay sum rules and factorization are presented. 
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I. INTRODUCTION 



Heavy quark symmetry implies that in the mg —>■ oo hmit matrix elements of the 
weak currents between a B meson and an excited charmed meson vanish at zero recoil (where 
in the rest frame of the B the final state charmed meson is also at rest). However, in some 
cases at order Aqcb/^q these matrix elements are not zero 0. Since most of the phase 
space for semileptonic B decay to excited charmed mesons is near zero recoil, Aqcd/'"^q 
corrections can be very important. This paper is concerned with rates for B semileptonic 
decay to excited charmed mesons, including the effects of Aqcd/''^q corrections. 

The use of heavy quark symmetry resulted in a dramatic improvement in our under- 
standing of the spectroscopy and weak decays of hadrons containing a single heavy quark, 
Q. In the limit where the heavy quark mass goes to infinity, mg — »• oo, such hadrons are 
classified not only by their total spin J, but also by the spin of their light degrees of freedom 
(i.e., light quarks and gluons), [^. In this limit hadrons containing a single heavy quark 
come in degenerate doublets with total spin, J± = s^i |, coming from combining the spin of 
the light degrees of freedom with the spin of the heavy quark, sg = |. (An exception occurs 
for baryons with se = 0, where there is only a single state with J = The ground state 
mesons with Q q flavor quantum numbers contain light degrees of freedom with spin-parity 

~ k ' giving a doublet containing a spin zero and spin one meson. For Q = c these 
mesons are the D and D*, while Q = b gives the B and B* mesons. 

Excited charmed mesons with sj'' = |^ have been observed. These are the Di and D2 
mesons with spin one and two, respectively. (There is also evidence for the analogous Q = b 
heavy meson doublet.) For q = u,d, the Di and D2 mesons have been observed to decay to 
TT and are narrow with widths around 20 MeV. (The Dgi and D*2 strange mesons decay 
to D^*^K.) In the nonrelativistic constituent quark model these states correspond to L = 1 
orbital excitations. Combining the unit of orbital angular momentum with the spin of the 
light antiquark leads to states with s'^" = and The I"*" doublet, (Dq,D*), has not 
been observed. Presumably this is because these states are much broader than those with 
sj^ = 1'''. A vast discrepancy in widths is expected since the members of the I"*" doublet 

of charmed mesons decay to D^*^ vr in an S'-wave while the members of the |^ doublet of 
charmed mesons decay to D^*^ tt in a D-wave. (An S'-wave Di D* n amplitude is allowed 
by total angular momentum conservation, but forbidden in the rriQ 00 limit by heavy 
quark spin symmetry 0.) 

The heavy quark effective theory (HQET) is the limit of QCD where the heavy quark 
mass goes to infinity with its four velocity, v, fixed. The heavy quark field in QCD, Q, is 
related to its counterpart in HQET, h^^\ by 



Q{x) = e" 



■imQV-x 



1 + -^ + ... 

2mn 



h^^^ , (1.1) 



where "^h^^ = h^^^ and the ellipses denote terms suppressed by further powers of Aqcd/^^q- 
Putting Eq. into the part of the QCD Lagrangian involving the heavy quark field, 

C = Q{ip- mo) Q, gives 

C = CnQET + SC + ... . (1.2) 
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The HQET Lagrangian || 

CnQ^^ = h^^hvDh^^'^ (1.3) 

is independent of the mass of the heavy quark and its spin, and so for Nq heavy quarks with 
the same four velocity v there is a U{2Nq) spin-flavor symmetry. This symmetry is broken 
by the order Aqcd/''^q terms in 6C, 



6C 



'2mQ 
where 



QiQ) , q(Q) 1 (1 



OZ = h'^' h^^^ , OL«i, = h^^^ I . (1.5) 
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The first term in Eq. ( |1.4D is the heavy quark kinetic energy. It breaks the flavor symmetry 
but leaves the spin symmetry intact. The second is the chromomagnetic term, which breaks 
both the spin and flavor symmetries. (In the rest frame, it is of the form JIq ■ -Bcoior, where 
JIq is the heavy quark color magnetic moment.) 

The hadron masses give important information on some HQET matrix elements. The 
mass formula for a spin symmetry doublet of hadrons H± with total spin J± = ± | is 

where the ellipsis denote terms suppressed by more powers of Aqcd/'^q and n± = 2J± + 1 
is the number of spin states in the hadron H±. The parameter A is the energy of the light 
degrees of freedom in the rriQ oo limit, Ai determines the heavy quark kinetic energyQ 

Af = {H^iv)\ h^^^ h(^^ \H^{v)) , (1.7) 

and A2 determines the chromomagnetic energy 

A? = {H^{v)\ I a^,G-^ h^^) \H^{v)) . (1.8) 

A and Ai are independent of the heavy quark mass, while A2 has a weak logarithmic depen- 
dence on vhq. Of course they depend on the particular spin symmetry doublet to which if± 
belong. In this paper, we consider heavy mesons in the ground state = | doublet and 
the excited s'^' = I"*" and doublets. We reserve the notation A, Ai, A2 for the ground 

state multiplet and use A', A'^^, A2 and A*, A*, A2 for the excited = |^ and |^ doublets, 
respectively. 



^Hadron states labeled by their four- velocity, v = pn/fnii, satisfy the standard covariant normal- 
ization {H{p'j,) \ H{ph)) = {2T,f 2Eh S^'ip'n - Ph)- 
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Particles 




m (GeV) 


1- 

2 


D, D* 


0-, 1- 


1.971 


1 + 
2 


Dl Dl 


0+, 1+ 


~ 2.40 


3 + 
2 


Dl, Dl 


1+, 2+ 


2.445 



TABLE I. Charmed meson spin multiplets {q = u,d). 



The average mass niH, weighted by the number of hehcity states 

_ n_mH_ + n+mH+ „s 

ruH = ■ (1.9) 

n+ + 71- 

is independent of A2. The spin average masses for the lowest lying charmed mesons is given in 
Table |. Identifying the 3^*^11 resonances observed at LEP with the bottom sj'' = I"*" meson 
doublet we can use their average mass, m'^ = 5.73 GeV 0, to determine the differences 
A' - A and A'^ - Ai: 

A' - A = s ill £1 ^ 0.39 GeV , 

rUb - rric 

A;-Ai = ^ B Bi £Zi ^ _o.23 GeV^ 1.10 

nib - TTLc 

The numerical values in Eq. ( |1.10| ) follow from the choices mi, = 4.8 GeV and ttIc = 1.4 GeV. 
To the order we are working, nib and in Eq. ( |1.10D can be replaced by and niD- 



This changes the value of A' — A only slightly, but has a significant impact on the value of 
A'^ — Ai. The value of A' — A given in Eq. (|1.1CI|) has considerable uncertainty because the 
experimental error on is large, and because it is not clear that the peak of the B^*^ti 
mass distribution corresponds to the narrow |^ doublet .| 

At the present time, A and Ai are not well determined. A fit to the electron energy 
spectrum in semileptonic B decay gives A ~ 0.4 GeV and Ai ~ — 0.2GeV^, but the 
uncertainties are quite large (A linear combination of A and Ai is better determined 
than the individual values.) 

The measured D* — D mass difference (142 MeV) and the measured — Di mass differ- 
ence (37 MeV) fix A2 = 0.10 GeV^ and A^ = 0.013 GeV^ Note that the matrix element of the 



^The Bsi and B*2 masses could also be used to determine A' — A from the relation 

A' - A = a; - A + (m'o - rn'o^ ) + 0(Aqcd "i./mj , 

with the analog of Eq. ( |1.1CI| ) used to fix A'^, — A, and m'^ — m'^^^ = — 114MeV. The B* has not been 
observed, but its mass can be determined from {tub* — rriBj — {tub* — mB) = {mc/mb) [{niD* — 
f^Ds) ~ {mn* — m-z))]. However, because of uncertainties in the Bsi and B*2 masses and the 
unknown order (Aqcb nT's/^c) term, this relation does not give a more reliable determination of 
A' - A than Eq. ( [rTO[ ). 
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chromomagnetic operator is substantially smaller in the excited sj'^ = | multiplet than in 
the ground state multiplet. This is consistent with expectations based on the nonrelativistic 
constituent quark model. In this phenomenological model, the splitting between members 
of a Q g meson spin symmetry doublet arises mostly from matrix elements of the operator 
sq ■ Sq6^{f), and these vanish for Qq mesons with orbital angular momentum. 

Semileptonic B meson decays have been studied extensively. The semileptonic decays 
B Deue and B D*eue have branching ratios of (1.8 ± 0.4)% and (4.6 ± 0.3)%, 
respectively 0, and comprise about 60% of the semileptonic decays. The differential decay 
rates for these decays are determined by matrix elements of the weak 6 — c axial- vector and 
vector currents between the B meson and the recoiling D^*^ meson. These matrix elements 
are usually parameterized by a set of Lorentz scalar form factors and the differential decay 
rate is expressed in terms of these form factors. For comparison with the predictions of 
HQET, it is convenient to write the form factors as functions of the dot-product, w = v-v', of 
the four-velocity of the B meson, v, and that of the recoiling D^*^ meson, v'. In the uiq — >■ oo 
limit, heavy quark spin symmetry implies that the six form factors that parameterize the 
B ^ D and B ^ D* matrix elements of the b c axial-vector and vector currents 
can be written in terms of a single function of w |Q. Furthermore, heavy quark flavor 
symmetry implies that this function is normalized to unity at zero recoil, w = 1, where 
the D^*^ is at rest in the rest frame of the B [|I^,|r2,|I| . The functions of w that occur 
in predictions for weak decay form factors based on HQET are usually called Isgur-Wise 
functions. There are perturbative as{mQ) and nonperturbative Aqcd/''^q corrections to the 
predictions of the tjiq oo limit for the B ^ D^*^ e semileptonic decay form factors. The 
perturbative QCD corrections do not cause any loss of predictive power. They involve the 
same Isgur-Wise function that occurs in the mq — > oo limit. At order Aqcd/''^q several new 



Isgur-Wise functions occur; however, at zero recoil, there are no Aqcb /rriQ corrections [112 



Expectations for the B ^ D^*^ eue differential decay rate based on HQET are in agreement 
with experiment [|1^]. 

Recently, semileptonic B decay to an excited heavy meson has been observed ||T^-|TB[. 
With some assumptions, CLEO [|1^] and ALEPH [|15| find respectively the branching ratios 
B{B^ Dieue) = (0.49 ±0.14)% and i3(5 ^ Die z>e) = (0.74±0.16)%, as well as the limits 
B{B D*eue) < 1% and B{B D*eue) < 0.2%. In the future it should be possible 
to get detailed experimental information on the B ^ DicVe and B D2 e differential 
decay rates. 

In this paper we study the predictions of HQET for B semileptonic decay to excited 
charmed mesons. This paper elaborates on the work in Ref. and contains some new 
results. In the infinite mass limit the matrix elements of the weak axial-vector and vector 
current between the B meson and any excited charmed meson vanish at zero recoil by heavy 
quark symmetry. Corrections to the infinite mass limit of order Aqcb/i^q and order asimq) 
are discussed. The corrections of order Aqcd/''^q are very important, particularly near zero 
recoil. 

Section II discusses the differential decay rate d^F/dw dcos6' for B {Di, D2) e z/e, where 
9 is the angle between the the charged lepton and the charmed meson in the rest frame of 
the virtual W boson. Corrections of order Aqcd /mq are included. At order Aqcd/^^q 
the B ^ Di zero recoil matrix element does not vanish and is expressible in terms of the 
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leading mq oo Isgur-Wise function, r, and A' — A (which is known in terms of hadron 
mass sphttings from Eq. ( 1.101 )). Away from zero recoil new Isgur-Wise functions occur, 



which are unknown. These introduce a significant uncertainty. The Aqcb/^q corrections 
enhance considerably the B semileptonic decay rate to the Di state, and for zero helicity the 
slope of dT{B Di eh'e)/dw at w = 1 is predicted. These corrections also reduce the ratio 
R = B{B D2eUe)/B{B Dieve) compared to its value in the irtq —>■ oo limit. The 
value of r at zero recoil is not fixed by heavy quark symmetry, and must be determined from 
experiment. The measured B — > DieUg branching ratio is used to determine (with some 
model dependent assumptions) |t(1)| = 0.71. The effects of perturbative QCD corrections 
are also discussed, with further details given in Appendix A. 

It is interesting to understand the composition of the inclusive B semileptonic decay rate 
in terms of exclusive final states. In Section III, the HQET predictions for the differential 
decay rates for B ^ D^eUe and B ^ D^eUe are investigated. The situation for the excited 

= 1^ multiplet is similar to the sj* = |^ multiplet discussed in Section II. Using a 
quark model relation between the leading rriQ oo Isgur-Wise functions for B decays to 
the sj^ = I"*" and sj' = charmed mesons (and some other model dependent assumptions), 
the rates for B Dq e z/g and B ^ Dleve are predicted. 

Section IV discusses the contribution of other excited charmed mesons to the matrix 
elements of the vector and axial- vector current at zero recoil. Only excited charmed hadrons 
with sY = \ , I and l"*^ can contribute. The I"*" and I"*" doublets are discussed in 
Sections II and III. This section deals with the \ and | cases, where the Aqcd/'^q 
corrections to the states from 5C give rise to non-vanishing zero recoil matrix elements. 

Section V examines other applications of our results. Using factorization, predictions 
are made for nonleptonic B decay widths to D2 vr, Di vr and to ir, Dq it. The importance 
of our results for B decay sum rules is discussed. Including the excited states dramatically 
strengthens the Bjorken lower bound on the slope of the B D^*^ eue Isgur-Wise function. 

Concluding remarks and a summary of our most significant predictions are given in 
Section VI. 



II. B ^Dieve AND B ^ D^eUe DECAYS 

The matrix elements of the vector and axial-vector currents {V^ = c'-f'^b and 
57^^75 &) between B mesons and Di or D2 mesons can be parameterized as 



{D,{v',e)\V^\Biv)) 



^rriD^ rriB 

^niD^ rriB 
{Dl{v',t)\A^^\B{v)) 



M e*^ + ifvX + fv.v'n (e* ■ v) 



7 ' 



y/niD* TUB 

{Dl{v',e)\V^^\B{v)) 
yjniD* rriB 



ikye^^^^^elXv^v'^, (2.1) 
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where the form factors and ki are dimensionless functions of w. At zero recoil {v = v') 
only the fy^ form factor can contribute, since v' dotted into the polarization (e*'^ or e*'^") 
vanishes. 

The differential decay rates can be written in terms of the form factors in Eq. (|2.1| ) . It is 
useful to separate the contributions to the different helicities of the Di and D2 mesons, since 
the Kqcd/itiq corrections effect these differently, and the decay rates into different helicity 
states will probably be measurable. We define 6 as the angle between the charged lepton and 
the charmed meson in the rest frame of the virtual W boson, i.e., in the center of momentum 
frame of the lepton pair. The different helicity amplitudes yield different distributions in 9. 
In terms oi w = v ■ v' and 9, the double differential decay rates are 



dw dcos 9 



d^r 



dw dcos 9 



+ (1 - 2r^w + rj) (1 + cos^ 9) [f^^ + {w' - - A cos 9 V^i^^ fy, /a 



1)3/2^1 sin^^ 



l~2r2W + rl) [(1 + cos2^) [k\^ + {w^ - 1)A;^] -4cos^ -IkA, h 



(2.2) 



where Fq = | m^/(1927r^), ri = niDi/mB, T2 = mD*/mB. The semileptonic B decay 
rate into any J 7^ 1 state involves an extra factor of w"^ — 1. The sin^ 9 term is the helicity 
zero rate, while the 1 + cos^ 9 and cos 6* terms determine the helicity A = ±1 rates. Since 
the weak current is V — A in the standard model, B mesons can only decay into the helicity 
|A| = 0, 1 components of any excited charmed mesons. The decay rate for |A| = 1 vanishes 
at maximal recoil, Wmax = (1 + '"^)/(2t'), as implied by the l — 2rw + r'^ factors above (r = ri 
or r2). From Eq. (|2.2| ) it is straightforward to obtain the double differential rate d^T/dwdy 
using the relation 



y = l 



rw — r 



Vw'^ — 1 cos 9 , 



(2.3) 



where y = 2Ee/mB is the rescaled lepton energy. 

The form factors fi and ki can be parameterized by a set of Isgur-Wise functions at each 
order in Aqcb /i^q- It is simplest to calculate the matrix elements in Eq. (|2.1| ) using the 



trace formalism |17,18|. The fields Py and P*'^ that destroy members of the sj^ 
with four-velocity v are in the 4x4 matrix 



I doublet 



Hv 



(2.4) 



while for sj^' = |^ the fields and P*^" are in 



The matrices H and F satisfy the properties 'f)Hy = = —Hy'^, = 
Fi:^, = Fi}v, = Q. 



(2.5) 
-F^l^, and 
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To leading order in Kqc-q/itiq and a^, matrix elements of the b ^ c flavor changing 
current between the states destroyed by the fields in Hy and F", are 

crb = h^^^r h^^^ = t{w) Tr {v^f;;, r h,} . (2.6) 

Here t{w) is a dimensionless function, and h^^^ is the heavy quark field in the effective theory 
(r is times the function T3/2 of Ref. [l^). This matrix element vanishes at zero recoil 
for any Dirac structure F and for any value of r(l), since the B meson and the {Di,D2) 
mesons are in different heavy quark spin symmetry multiplets, and the current at zero recoil 
is related to the conserved charges of heavy quark spin-flavor symmetry. Eq. ( p.6|) leads to 
the rriQ —>■ 00 predictions for the form factors and ki given in Ref. ||19|| . 

At order Aqqd /mq, there are corrections originating from the matching of the b ^ c 
flavor changing current onto the effective theory, and from order Aqcd/'"^q corrections to 
the effective Lagrangian. The current corrections modify the first equality in Eq. (p.6|) to 

For matrix elements between the states destroyed by the fields in F^, and H^j, the new order 
Aqcd/'^q operators in Eq. (|2.7| ) are 

hif iDx V /if = Tr [S^l F H,] , 

h!^) V i^^ h^^ = Tr {^S TH,]. (2.8) 
The most general form for these quantities is 



(Q) 



(Q) , (Q) / , (Q) 



9aX ■ (2.9) 



The functions depend on w, and have mass dimension one.0 They are not all independent. 
The equation of motion for the heavy quarks, (f ■ D) h[^^ = 0, implies 

rf)+ti;rf -rf +rf =0. (2.10) 

Four more relations can be derived using 

td, {h^) F /il")) = {kv, - k'v',) h^) F /il^) , (2.11) 

which is valid between the states destroyed by the fields in F^, and H^. This relation follows 
from translation invariance and the definition of the heavy quark fields It implies that 



^Order Aqcd/^t-c corrections were also analyzed in Ref. [20|. We find that r4 (denoted ^4 in 
does contribute in Eq. ( |2.8| ) for T = 7a F, and corrections to the Lagrangian are parameterized by 
more functions than in 
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c(c) I qib) 



T . 



(2.12) 



Eq. ( p^.l2D gives the following relations^ 



-ic) 



-ib) 



At, 
-A'r, 

0, 
0. 



(2.13) 



These relations express the 'T-^'^'s in terms of the Tj'^'''s. Furthermore, combining Eqs. (p.lO|) 
with (|27[3| ) yields 
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(c) , (c) 



(W 



l)(r- 



ri"^) - {wA' -A)t. 



(2.14) 



All order Aqcd/'"^q corrections to the form factors coming from the matching of the QCD 
currents onto those in the effective theory are expressible in terms of A r and A' r and two 



functions, which we take to be and t^' . From Eqs. (|2.8| ) and (2^) it is evident that 
only rf^'' can contribute at zero recoil. Eq. (|2.14|) determines this contribution in terms of 
r(l) and measurable mass splittings given in Eq. 



LIOD 



-(b). 



-(c) 



(1) = (A'-A)r(l) 



(2.15) 



Note that with our methods Eq. (|2.15| ) cannot be derived working exclusively at zero 
recoil. At that kinematic point, matrix elements of the operator h^f> F /i^^^ vanish between a 



B meson and an excited charmed meson, and so Eq. ( p.ll|) only implies that rf' + rj"'' = 0. 



Sb) 



Eq. ( p.l5| ) relies on the assumption that the Tj^\w) are continuous at w = 1. 



Next consider the terms originating from order Aqcd/''^q corrections to the HQET 
Lagrangian, 5C in Eq. ( |1.4| ). These corrections modify the heavy meson states compared to 
their infinite heavy quark mass limit. For example, they cause the mixing of the Di with 



the = 1"*" member of the s 



1 + 



doublet. (This is a very small effect, since the Di 



is not any broader than the D^O For matrix elements between the states destroyed by the 
fields in F", and the time ordered products of the kinetic energy term in 5C with the 
leading order currents are 

z / d^xT{0S,,,(x) [/.J)f/.W](0)} = r^£?Tr{t;.F;Ff/.}, 

T {0£i,,(x) [/.J) F h^)] (0) } = r/(? Tr {v^F^^, F if.} . (2.16) 

These corrections do not violate spin symmetry, so their contributions enter the same way 
as the mg 00 Isgur-Wise function, r. 



^In Ref. p two out of these four relations were obtained (only those two were needed to get 
Eq. ( |2.15| )). We thank M. Neubert for pointing out that there are two additional constraints. 
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For matrix elements between the states destroyed by the fields in and H^, the time 
ordered products of the chromomagnetic term in 6C with the leading order currents are 



d^x 



;T{0 



mag,!)' 



(b) ( 

mag,i) v 



Tr \ n% F:, la 



(2.17) 



The most general parameterizations of T^^*^) are 



n. 



aafi 



(c) , (c) , (c) 



Tlib) (0) , (0) / , (0) I 



(2.18) 



Only the part of 7^^*^^ antisymmetric in a and (3 contributes when inserted into Eq. ( ^.17^ . 
The functions rji depend on w, and have mass dimension one. Note that g^alp is dependent 
on the tensor structures included in Eq. (|2.18| ) for matrix elements between these states. 
For example, for the Aqcd / corrections the following trace identity holds 



Tr 







(2.19) 



All contributions arising from the time ordered products in Eq. ( p.l7| ) vanish at zero recoil, 
since v^F" = and ^^(1 + ^)(j"'^(l + ^) = 0. Thus we find that at zero recoil the only 
Aqcb/^q corrections that contribute are determined by measured meson mass splittings 
and the value of the leading order Isgur-Wise function at zero recoil. 

The form factors in Eq. ( |2.1| ) depend on rjl''^ only through the linear combination r/^ = 
rjl^J + QVi'^ — 2(w — 1) 7^2^^ + . Denoting Eq = l/{2mQ) and dropping the superscript on 
Tj*''^^ and ini\^\ the B —>■ DieUe form factors are [@ 

VefA = -{w + 1)t - Ebiiw - 1)[(A' + A)r - {2w + l)n - ra] + {w + l)rjb} 
-e,[A{wA' - A)r - 3{w - l){n - r^) + {w + l)(r/ke - 2r]i - Sr/g)] , 
(1 - w^)t - Ebiw'' - 1)[(A' + A)r - {2w + l)ri - T2 + Vb] 
-e,[4{w + 1)(«;A' - A)r - {w^ - l){3n - 3r2 - r]^, + 2r]^ + Sr/g)] , 
-3r - 3£5[(A' + A)r - {2w + l)n - ra + Vb] 
-ec[{4:W - l)ri + 5r2 + Sr^kc + lOr/i + A{w - l)r]2 - Sr^a] , 
{w - 2)t + £h{{2 + w)[{A' + A)r - {2w + l)n - ^2] - (2 - w)7]i,} 
+e,[4:{wA' - A)r + (2 + w)n + (2 + 3w)t2 

+ {w- 2)r7kc - 2(6 + w)r]i - A{w - l)r]2 - {3w - 2)773] • (2.20) 

The analogous formulae for B D2 e z/g are 

kv = -T- EbliA' + A)r - {2w + 1)ti -T2 + 775] - £c{ti -T2 + rjy^^ - 2r]i + 773), 
kA, = -(1 + w)t - 6b{{w - 1)[(A' + A)r - {2w + l)ri - T2] + (1 + w)r]b} 

-Ecliw - l)(ri - T2) + {W+ l)(77ke - 2r]i + 773)], 

kA, = -2e, in +7]2), (2.21) 
kA3 = r + 6b[{A' + A)t - {2w + l)ri - T2 + r]b] - edn + T2 - rj^e + 2771 - 2r]2 - 773). 
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Recall that fy^ determines the zero recoil matrix elements of the weak currents. From 
Eqs. ([2:201 ) it follows that 



\/6 fv, (1) = -8 £e (A' - A) r(l) . (2.22) 

The allowed kinematic range for B — >■ DieUe decay is 1 < if < 1.32, while for B — *■ 
D2 e Vf. decay it is 1 < w < 1.31. Since these ranges are fairly small, and at zero recoil there 
are some constraints on the Aqcd/'^q corrections, it is useful to consider the decay rates 



given in Eq. (2^) expanded in powers of w — 1. The general structure of the expansion of 
dr/dw is elucidated schematically below. 



dr 



(A=0) 



dw 



u;- l)°(0 + 0£ + £2 + e3 + .. 

+ (u; - 1)^ (0 + £ + + . . . ) + (ly - 1)2 (^1 + £ + . . . ) + 



dr 



(|A|=1) 



dw 



jr(|A|=o,i) 



w 



+ Oe + £^ + e^ + ... 



+ {w-iy {1 + £ + ...) +{w- if (l + e + ... 



dw 



[w 



\3/2 



l)°(l+£ + ...) + («;-l)i(l+£ + . ..) + ... 



(2.23) 



Here denotes a term of order (Aqcd/'^q)"- The zeros in Eq. ( p.23| ) are consequences 
of heavy quark symmetry, as the leading contribution to the matrix elements of the weak 
currents at zero recoil is of order Aqcd/^^q- Thus, the Di decay rate at w = 1 starts out at 
order Aqqy)/itlq. Similarly, from Eq. (1^) it is evident that the vanishing of /vi(l) in the 

niQ 00 limit implies that at order w — 1 the D[^~^^ rate starts out at order Aqcd/tuq. 
The D2 decay rate is suppressed by an additional power of — 1, so there is no further 
restriction on its structure. 

In this paper we present predictions using two different approximations to the decay 
rates. In approximation A we treat w — 1 as order Aqqb/^q and expand the decay rates 
in these parameters. In approximation B the known order Aqcb/^q contributions to the 
form factors are kept, as well as the full w-dependence of the decay rates. 

Expanding the terms in the square brackets in Eq. (|2.2| ) in powers of w — 1 gives 



d^r 



d..dcos« r„r^(l).?v/^i: (--!)"{ -="4' 

+ (1 - 2riw + rl) [(1 + cos^ 6) t^"^ - 4 cos ^ 

+ (1 - 2r2W 



(2.24) 



d^r 



dw dcos 9 



w 



1)3/2 ^ 



W 



1)N- sin2^4") 



(1 + cos^ e) t^"^ - 4 cos ^ Vw^ - 1 



(n) 
U2 



(We do not expand the factors of \/w'^ — 1 that multiply cos 6*). The subscripts of the 
coefficients s, t, u denote the spin of the excited D meson, while the superscripts refer to the 
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order in the w — 1 expansion. The terms proportional to cos^ only affect the lepton 
spectrum, since they vanish when integrated over 6. 

Eqs. ( p.2D , ( |2.20|) , and ( |2.21| ) yield the following expressions for the coefficients in the Di 
decay rate in Eq. (|2.24|) , 



.(0) 

'1 

.(1) 



,(2) 



(0) 



.(2) 



32£2(l-ri)'(A'-A)' + ... , 
32sc (1 - r^) (A' - A) + . . . , 

8{l + nY + ... , 
32el (A' - A)2 + . . . , 

4 + See [4(A' - A) + r)ke - 2fii - 8173] + Sefer),, + . . . 

8 (1 + f') + . . . , 
Ssc (A' - A) + ... , 
2 + ... . 



(2.25) 



For the decay rate into D2 the first two terms in the w — 1 expansion are 



JO) 

AO) 



4(1- [I + 2eb fjb + 2e, {fj^, - 21)^ + fj^) 

4(l-r2)'(l + 2f') + ... , 

4 + S^b r)b + (57ke - 2171 + f/g) + . . . , 

2(3 + 4f') + ... , 

2 + ... . 



+ . . . 



(2.26) 



In Eqs. ( p.25| ) and (p.26| ) the functions r, r' = dr/dw, and r/j are all evaluated at w = 1, and 
the functions with a hat are normalized to r(l) (e.g., fji = ?7j/r(l), f = r'/r(l), etc.). The 
ellipses denote higher order terms in the Aqcd / ""^q expansion. The -u terms are suppressed 
by y/w"^ — 1 compared to s[^^ and t\"'\ therefore we displayed the m's to one lower order than 
the s and t coefficients. (Note that u^i^ also starts out at order Aqcd/^^q as a consequence 
of the vanishing of fvii^) in the rriQ 00 limit, as it was shown for s^i^ after Eq. ( p.23|) .) 

The order Aqcd / "^q terms proportional to A' — A are very significant for the Di decay 
rate. The decay rate into D2 does not receive a similarly large enhancement from order 
Aqcd/'^q terms proportional to A' — A. The coefficients s^^^ and t^^'' are independent of A' 
and A to the order displayed in Eq. (f 



The values of Si and t[ are known to order Aq^j-j/mQ, and Si and Ui are known 

to order Aqcd/'^q- At order Aqcd/^^q? the only unknowns in S2°\ 4°^ the fji 
functions that parameterize corrections to the HQET Lagrangian. The remaining coefficients 

^(2) ^(2) ^,(1) Jlj ^(1) ^,(0)^ 



in Eqs. and ( ^1261) (i.e., s\ ', t 



'1 ) ''1 ) ^1 ) 



are known in the infinite mass 



limit in terms of f'(l), the slope of the niQ — > cxo Isgur-Wise function at zero recoil. At order 
Aqcd/'"^q, these six coefficients depend on the unknown subleading and r]i functions. 

The values of r', iqf^^ and ri 2 that occur in Eqs. (|2.25|) and ( p. 261) are not known (r, only 
appears in the terms replaced by ellipses). 771^^3, which parameterize time ordered products 
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of the chromomagnetic operator, are expected to be small (compared to Aqcd), and we 
neglect them hereafter. This is supported by the very small D2 — Di mass splitting, and 
the fact that model calculations indicate that the analogous functions parameterizing time 
ordered products of the chromomagnetic operator for B D'^*'> ev^ decays are small . 
On the other hand, there is no reason to expect ri 2 and t]^"* to be much smaller than about 
500 MeV. Note that the large value for A'^^ is probably a consequence of the Di and D2 being 
P-waves in the quark model, and does not necessarily imply that significantly distorts 
the overlap of wave functions that yield //^^^ 

Even though ec{^' — A) ~ 0.14 is quite small, the order Aqcd/"^q correction to 
proportional to ec(A' — A) is as large as the leading mq — > 00 contribution. This occurs 
because it has an anomalously large coefficient and does not necessarily mean that the 
Aqcd/'^q expansion has broken down. For example, the part of the Aq^j-j/m^ corrections 

that involve A', A, and r'(l) affect s^^^ by (21 + 10f')%, and t^i^ by (44 + 15f')% (using 
A = 0.4 GeV 0). These corrections follow from Eq. (|2.20|) , but they are neglected in 
Eq. (|2.25| ) (i.e., approximation A), because there are other order Aqcp/mQ effects we have 
not calculated. 

As the kinetic energy operator does not violate spin symmetry, effects of t]^^ can be 
absorbed into r by the replacement of r by f = r + V\e + ^fe V\I ■ This replacement 
introduces an error of order Aqcj-j/mg, in t^^^ , etc. But due to the presence of large Aqcd / 
corrections, the resulting ^^q^/^'q error is also sizable, and is expected to be more like an 
order Aqcd/'^q correction. Hereafter, unless explicitly stated otherwise, it is understood 
that the replacement r — >■ r is made. But we shall examine the sensitivity of our results 
to r^ke (assuming it has the same shape as r). 

In approximation A we treat w — 1 as order Aqcd/'^q 0! and keep terms up to order 
(AQCD/mQ)2— in 4"^ and t^"^ (n = 0, 1, 2) in Eq. Q^ ), and up to order (Aqcd/^q)'"" in 
s^"^ and t^2^ (n = 0, 1) in Eq. ([2:261) . Since the uf'' are suppressed by — 1 compared 
to s|"'' and we keep Uj-"'' to one lower order than the s and t coefficients, i.e., to order 
(AqcdMq)^"'' {n = 0, 1) for B ^ Di decay and order (Aqcd/?t^q)" {n = 0) for B ^ DI 
decay. The terms included in approximation A are precisely the ones explicitly shown in 
Eqs. (|2.25|) and ( p. 261) . This power counting has the advantage that the unknown functions, 
Ti and r2, do not enter the predictions.^ Neglecting higher order terms in the w — 1 expansion 
in approximation A gives rise to a sizable error for the B ^ DicVf. decay^. The order {w—lY 
term is important for the decay into helicity zero Di in the itlq 00 limit, since the helicity 
zero rate (which, as we shall see, dominates over the helicity one rate) starts out at order 
(w - 1)2 as shown in Eq. (^^). 

In approximation B we do not expand the decay rates in powers of w — 1. We keep the 



^Approximation A differs from our discussion in Ref. Q only in the separation of the different 
lielicity states of the excited charmed mesons, and keeping the 1 — 2rw + factors for the helicity 
one states as well as the {w'^ — 1)^/^ terms for the D2 rates unexpanded. 

^We thank A. Le Yaouanc for pointing out the importance of these terms. 
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Aqcd / "^Q corrections to the form factors that involve A' and A and examine the sensitivity 
of our results to the corrections involving ti and T2 (assuming that they have the same shape 
as r, which is not a strong assumption). This approximation retains some order K^Qj^/rriQ 
terms away from zero recoil in the differential decay rates. Furthermore, a linear form for the 
Isgur-Wise function is assumed, t{w) = r(l) [l-\-f'{w — l)]. The uncertainty in the Aqcd/'^q 
corrections is parameterized by the functions ti^2{w). A different choice of ti^2{w) changes 
what is retained by terms involving K/mq and A'/mg. In an approximation, which we shall 
refer to as Bi, we set ri = r2 = in Eqs. ( |2.2CI|) and (|2.21|) . (This is identical to saturating 
the first two relations in Eq. ( p.l3|) by rf^, i.e., setting r^''^ = Ar and T^f* = — A'r.) An 
equally reasonable approximation, which we refer to as B2, is given by setting ti = At and 
T2 = -K't in Eqs. (|2]20|) and {^2^. (This is identical to setting r}^ = 0.) If the first 
two relations in Eq. ( ^.131 ) are taken as hints to the signs of ti and T2, then the difference 
between approximations Bi and B2 gives a rough estimate of the uncertainty related to the 
unknown Aqcd/^^q corrections. When our predictions are sensitive to Ti and T2, we shall 
vary these in a range larger than that spanned by approximations Bi and B2. Note that 
the infinite mass limits of Bi and B2 coincide. Predictions of approximation A are within 
the spread of the approximation B results, except for those that depend on the hehcity zero 
Di rate. In that case, including the order {w — 1)^ term in the infinite mass limit alone, 
s^i^ = 8 (l + ri)^ (l + 2f'), would bring the approximation A results close to approximation B. 

Eqs. ( ^.25| ) and (|2.2(j| ) show that the heavy quark expansion for B decays into excited 
charmed mesons is controlled by the excitation energies of the hadrons. A' and A. For highly 
excited mesons that have A' comparable to rric, the l/mq expansion is not useful. For the 

~ 2~^ doublet €cA' ~ 0.3. However, near zero recoil only Ec (A' — A) ~ 0.14 occurs at 
order AqcdMq- 

The expressions for the decay rates in terms of form factors in Eq. ( p.2|) imply that 
one form factor dominates each decay rate near zero recoil, independent of the helicity 
of the Di or D2 {fvi for Di and fc^i for D2). Thus, to all orders in the Aqcd/'^q 
expansion, Si^^ /ti'^ = (1 — ri)^, and S2°V^2°'' = (1 ~ ^2^- This implies that for 



B Di decay lim^ 



[(drgr°Vd«;)/(drK'-^VduO 

lim^^i (dr^r^Vdw) I {dT'j^}~^'' /dw) =2/3. Note that the first of these ratios would van- 
ish if the rates were calculated in the mq — > 00 limit. In that case fv^iX) = 0, so the ratio 
of helicity zero and helicity one B ^ Di rates is determined by the other form factors at 
zero recoil. 



1/2, and for B ^ D2 decay 



A. Predictions 

The relationship between sf*^ and s^^^ implies a model independent prediction for the 
slope parameter of semileptonic B decay into helicity zero Di. This holds independent of 
the sub leading Isgur-Wise functions that arise at order Aqcd/'^q- The semileptonic decay 
rate to a helicity zero Di meson is 



dr 



(A=0) 
Di 



128 



dw 



Fo rl (1 - n f V^I^^r\l) el (A' - A)^ [l - p^^ {w - 1) + 



(2.27) 
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Approximation 






„(A=0) /t^ /i\ 
^Dl I^Bl t(1) 


■ 6.0 X 10-3 ■ 


1/2 


B{B Dieue). 






0.93 


0.88 


0.64 


0.92 


Boo 


1.65 


0.80 


0.66 


1.24 


A 


0.40 


0.81 


0.64 


0.60 


Bi 


0.52 


0.72 


0.63 


0.71 


B2 


0.67 


0.77 


0.64 


0.75 



TABLE II. Predictions for various ratios oi B ^ Di e Ue and B 
described in the text. The extracted value of r(l) is also shown. Aqo and Bqo denote the tuq 
limits of approximations A and B. These results correspond to f = r'(l)/r(l) = —1.5. 



D2 e Ve decay rates, as 



00 



where the slope parameter p|,^ for hehcity zero Di has the value 



1 + Ti 2mc 



1 - ri A' - A 



+ 0{l). 



(2.28) 



Since the decay rate at zero recoil is suppressed, p|,^ is of order mg/AqcD- Note that this 
slope parameter is negative. 

Recently the ALEPH |T3[ and CLEO [|TB| Collaborations measured, with some assump- 
tions, the B —>■ DieUe branching ratio. The average of their results is 

B{B ^ Die Ue) = (6.0 ± 1.1) x 10~^ . (2.29) 

The B D2 e branching ratio has not yet been measured; CLEO set the limit 
B{B ^ D*e Ue) < 1% PI, while ALEPH found B{B D*e Ve) < 0.2% p . 

Predictions for various quantities of experimental interest are made in Table || using 
A' - A = 0.39 GeV, A = 0.4 GeV, tb = 1.6 ps, \Va\ = 0.04, = 1.4 GeV, = 4.8 GeV. 
Keeping mj, — rric fixed and varying rric by ±0.1 GeV only affects our results at the few 
percent level. These predictions depend on the shape of the Isgur-Wise function. In our 
approximations this enters through the slope parameter, f = r'(l)/r(l), which is expected 
to be of order —1. We shall quote results for the "central value" f = —1.5, motivated by 
model predictions p^-p5| , and discuss the sensitivity to this assumption. For B ^ DieUe 
decay we use ri = 0.459 and 1 < w < 1.319, whereas for B ^ D^eUe decay r2 = 0.466 and 
l<w < 1.306. 

The order Aqcb/itT'Q corrections are important for predicting 



R 



B{B ^D*eue 
B{B ^ Dieue 



(2.30) 



In the niQ —>■ 00 limit R ~ 1.65 for f' = —1.5 (this is the B^o result in Table H). The sizable 
difference between approximations A and B is mainly due to the order {w~l)^ contribution to 
the helicity zero Di rate. For f = —1.5 this term by itself would shift the approximation A 
result for R from 0.40 to 0.49 and the A^o prediction from 0.93 to 1.65. The Aqcd/"^q 
correction to the form factors yield a large suppression of R as shown in Table H and Fig. |l|a. 
Fig. |l|a also shows that R is fairly insensitive to f'. The difference of the Bi and B2 results in 
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FIG. 1. Fig. la shows R = B{B D^ei?e)/B{B Dieue) as a function of f'. The dotted 
curve is the niQ — > oo hmit (Bqo), sohd curve is approximation Bi, dashed curve is B2. Fig. ^3 
shows i? as a function of ti(= ti/t) for r2 = (sohd curve), and as a function of T2 for f 1 = 
(dashed curve). Note that the scales in Fig. ||a and |l|b are different. 



Table |T|and Fig. |T]a shows that R is sensitive to the unknown Aqcb/itT'Q corrections, ti and 
T2. In Fig. [I|b we plot R in approximation B as a function of fi setting 7^2 = (solid curve), 
and as a function of f2 setting fi = (dashed curve). Fig. |l]b shows that R is fairly insensitive 
to r2, whereas it depends sensitively on ri. In the range —0.75 GeV < fi < 0.75 GeV, R 
goes over 0.27 < R < 1.03. This suppression of R compared to the infinite mass limit is 
supported by the experimental data. (It is possible that part of the reason for the strong 
ALEPH bound B{B ^ D*eiyeX)x B{D* ^ L)(*)7r) < (1.5 - 2.0) x 10"=^ |1| is a suppression 
of B{D* D'-*h) compared to B{Di D*n).) 

The prediction for the fraction of helicity zero Dis in semileptonic B ^ Di decay, 
F^~°^/F£i^, is surprisingly stable in the different approximations (see Table H). The weak 
dependence of this ratio on f is well described in approximation B for |1.5 + f'l < 1 
by adding 0.05(1.5 + f'). The dependence on ri is at the 0.01 level, while the T2- 
dependence is — 0.07f2/GeV. This is why the B2 result for this quantity is 0.05 larger 
than the Bi prediction. A linear dependence of (dF^~°''/dw) /(dF^^/dw) on w between 



lim„ 



(dFg=°Vd^)/(dF^Jdt.) 



(dFg=°Vd^)/(dF^,/d^) 



1 aX w 



1/3 and 

is consistent with our result. 

A similar prediction exists for the fraction of helicity zero -Dg's in semileptonic B 



D* 



decay. As can be seen from Table |T|, it is again quite stable. The dependence on f in 
approximation B is given by adding 0.04(1.5 + f'). However, Fl^r^VF^). is sensitive to both 
Ti and T2 at the (10 — 20)% level, and the small difference between the Bi and B2 predictions 
for this quantity in Table ^ is due to an accidental cancellation. The prediction for the w 



dependence of (dF^. /dw) j (dF£)*/dw) between lim^^i (dF^* jdvS) j (dF£)*/dw) =2/5 



and 



(dFgfVd^) 



.(A=0) 



'(dF^./d^) 



1 at w 



in this case is not linear. 
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FIG. 2. Fig. |2|a shows the extracted value of r(l) as a function of f' in approximations Bqo, 
Bi, and B2. The notation is the same as in Fig. |^. Fig. |2|b shows the dependence of t(1) on fi for 
f' = —1 (dashed curve), f' = —1.5 (sohd curve), f = —2 (dash-dotted curve). 



The predictions considered so far do not depend on the value of t(1), but r(l) affects 
some results that we discuss later. r(l) can be determined from the measured B —y Diev'e 
branching ratio using the expressions in Eq. ( |2.24D and ( [2.251 ). Using approximation Bi and 
f = —1.5, we obtain 



r(l) 



6.0 X 10- 



1/2 



0.71 



(2.31) 



The extracted value of r(l) is plotted in Fig. ||a in approximations Bqo, Bi, and B2 as 
functions of f. The suppression of r(l) compared to the infinite mass limit indicates that 
the order Kqcd/itiq corrections enhance the semileptonic B ^ Di width by about a factor 
of three. In approximation B the value of r(l) changes by less than 0.01 as T2 is varied 
in the range —0.75 GeV < 7^2 < 0.75 GeV, but r(l) is sensitive to ri at the 15% level. In 
Fig. Hb we plot r(l) as a function of fi for f' = — 1 (dashed curve), f' = —1.5 (solid curve), 
and f = —2 (dash-dotted curve). For ri > (such as approximation B2) r(l) is enhanced 
compared to the Bi value of 0.71. 

The value of r(l) in approximation B is larger than that in approximation A. Most of 
the difference arises from the inclusion of the order {w — 1)^ term, which reduces the 
theoretical expression for the helicity zero B DicVe rate (for f' < —0.5), resulting in 
an increase in the value of r(l) needed to accommodate the measured rate. For f = —1.5 
this term by itself would shift the approximation A result from 0.60 to 0.66, and the Aqo 
prediction from 0.92 to 1.22. The ISGW nonrelativistic constituent quark model predicts 
t(1) = 0.54, in rough agreement with Eq. (|2.31| ) p^JT9[] . (For some other quark model 
predictions, see, e.g., Ref. |]2^p5| . QCD sum rules can also be used to estimate r, see, e.g., 
Ref i2F 
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Approximation 



5(r,./r,,) «-7r,,) ^(rgfVr,,) SrH) 



6.0 X 10- 



_BiB Dieue) 



1/2 



-0.68 
-1.63 



0.10 
0.19 



0.02 
-0.003 



-0.26 
-0.32 



-0.22 
-0.55 
-0.68 



0.04 
0.06 
0.07 



0.05 
-0.02 
-0.05 



-0.24 
-0.32 
-0.33 



TABLE III. Order Og and as{A.QCB/iT^Q) corrections to the predictions in Table |I^ for 
= —1.5. These numbers should be multiplied by as{^/ri\jn^) / n to get the corrections to Table ||. 



The ALEPH and CLEO analyses that yield Eq. (1^291) assume that B DieUeX is 
dominated by -B — DieUe, and that Di decays only into D*7i. If the first assumption 
turns out to be false then r(l) will decrease, if the second assumption is false then r(l) will 
increase compared to Eq. ( p^.31| ). 

The predictions discussed above would change if we had not absorbed into r the time 
ordered product involving the kinetic energy operator. As discussed earlier (in the paragraph 
preceding the description of approximation A), the replacement of r by r = r+Ecfjl^J +eb'r]]^J 
introduces an error, which is formally of order AQQY^/mQ. Absorbing r^ke into r almost 
fully eliminates the ?7ke dependence of the D2 rate. For the Di rate, however, absorbing 
r^ke into r generates at order Aqc'd/''^q ^ formally suppressed but numerically sizable ?7ke 
dependence. This T^ke dependence is more like a typical Aqcd/'"^q correction, since the 
Aqcd/'"^q current corrections are as important as the infinite mass limit for the Di rate. 
Keeping 17^^^ = ij^^ /t explicit in the results, the total B —>■ Di semileptonic rate in units of 

ror2(l) is 0.033(1 + 1.1 + •••), while the B ^ D* rate is 0.017 (1 + 2.0 + •.•)• 

From these expressions it is evident that, for — 0.75GeV < r/ke < 0.75 GeV, r(l) changes 
only by ±15%, while R has a larger variation. In the future this uncertainty will be reduced 
if differential spectra can also be measured besides total rates in B ^ Di, D2 decays. Note 
that r/ke does not enter into predictions for the B DieUe decay rate near zero recoil. 

Order as corrections to the results of this section can be calculated in a straightforward 
way, using well-known methods. Details of this calculation are given in Appendix A. The 
order corrections to the results shown in Table y are given in Table |in[ These are smaller 
than the uncertainty in our results from higher order terms in the Aqcd / '"^q expansion that 
have been neglected. The corrections are most significant for R = r/j^^r^i^ and r(l) in 
approximation B; the central values of these quantities are reduced by about 9% and 4%, 
respectively. Some of these corrections depend sensitively on f, but they remain small 
for > f' > —2. For the remainder of this paper, we neglect the small corrections. 

Our predictions for the single differential B — > {Di, D2) eug spectra follow from 
Eqs. ( ^ ), (|;25D, and (^^. dT/dw is given by integrating Eqs. over dcos^. This 

amounts to the replacements sin^ 6* — > 4/3, (1 + cos^ 9) —> 8/3, and cos6' 0. Thus dT/dw 
is trivial to obtain using either approximations A or B. The electron energy spectra are ob- 
tained by expressing cos^ in terms of y (where y = 2Ef./mB is the rescaled electron energy) 
using Eq. (|2.3|), and integrating tu over [(1 — y)^ + r^]/[(2r(l — y)] <w < (l + r^)/(2r). They 
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depend on the coefficients u^f^ which did not enter our results so far. 

In Fig. ^the electron spectrum for ^ Z^i e i>e is plotted in units of Fq r^(l). Figs. ^ and 
are the spectra for helicity zero and helicity one Di, respectively. In these plots f' = —1.5. 
The dotted curve shows the tuq — >• oo limit (Bqo), the solid curve is approximation Bi, the 
dashed curve is B2. Note that the kinematic range for ?/isO<?/<l — r^. Near ?/ = and 
y = 1 — r'^ the spectrum is dominated by contributions from w near Wmax- In this case, we 
expect sizable uncertainties in our results, for example, from unknown terms that occur in 
the uf^ terms in Eq. ( p.25|) at a lower order than in the s and t coefficients. Fig. ^ shows the 
large enhancement of the Di rate due to order Aqcd/^q corrections, and that the difference 
between approximations Bi and B2 is small compared to this enhancement. In Figs. ^ and 
^ we plot the electron spectrum for B D2eUe for helicity zero and helicity one D2, 
respectively. In this case the Aqcb/^q corrections are less important. 



III. B ^ Dleve AND B ^ D^eue DECAYS 



The other low lying states above the D^*^ ground states occur in a doublet with sj'' 



1 + 



These states are expected to be broad since they can decay into D^*^ tt in an 5- wave, unlike 
the Di and D2 which can only decay in a D-wave. (An S-wave decay amplitude for the Di 
is forbidden by heavy quark spin symmetry 0.) This section repeats the analysis of the 
previous section for these states. Since the notation, methods, and results are similar to 
those used in Sec. II, the discussion here will be briefer. 

The matrix elements of the vector and axial currents between B mesons and Dq ot 
mesons can be parameterized by 



{D;{v')\V'^\Biv)) = 
{D*,{v')\A'^\B{v)) 



{Dl{y',t)\ \B{y)) 



{Dl{v',e)\A^\B{v)) 
^rriDi rriB 



gv,e*'^+{gvX + 9v,v'n {e* ■ v) 



(3.1) 



where gi are functions of w. At zero recoil the matrix elements are determined by ^'+(1) and 
gy^{l). In terms of these form factors the double differential decay rates for B — >■ D^eVf, 
and B ^ DlePe decays are 



d^Fn. 



dw dcos 9 

d^r^> 

dw dcos 9 



3Fo {w' - If/' sin^ ^(^i + r*,) g+ - {1 - r*) g 
3Fo rf ^f^F^ I sin^ 9 \{w - r\)gy, + {w" - \){gy, + r\gy, 
+ (1 - 2r> + rf) [(1 + cos^ 9) \g\ + {w" - \)g'A - ^^0^9 ^J^^I^gy, gA 



(3.2) 



where Fq = G\ | Vchp m^/(1927r^), Tq = rnoiImB and r\ = mDi/mB- 
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FIG. 3. Electron spectrum for i? ^ Di e in units of Tq t^(1) for f' = —1.5. Figs. ^ and ^3 
are the spectra for helicity zero and helicity one Di, respectively. Dotted curves show the mq oo 
limit (Boo), solid curves are approximation Bi, dashed curves are B2. 
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We follow the previous section to obtain expressions for the form factors gi in terms of 
Isgur-Wise functions to order Aqcd/'"^q- The fields Py and P*'^ that destroy members of 
the sj'' = doublet with four-velocity v are in the 4x4 matrix 

K = ^[p:^1,1, + P.]. (3.3) 

This matrix K satisfies -^K^ = Ky = Ky-^. In the infinite mass limit matrix elements of the 
leading order current operator are |T^ 

hif r /iW = ({w) Tr {k,> r H,} (3.4) 

Here C('"^) is the leading order Isgur-Wise function {( is twice the function ti/2 of Ref. [0). 
Since the (Dq,DI) states are in a different spin multiplet than the ground state, (?+(!) = 
(7vi(l) = in the infinite mass limit, independent of C(l)- 

The order Aqcd/'^q corrections to the current can be parameterized as 

iDx r /li") = Tr {5f K,, r if,} , 

r i^^ = Tr {5f K,, TH,]. (3.5) 
This is the analogue of Eq. ( p.8|) , except that in the present case 

The functions C,\^^ (w) have mass dimension one. The heavy quark equation of motion yield 

t^d^^+d'^^+cf = 0, 

d'^+^d'^-d'^ = o. (3.7) 

Eq. ( pID implies Si""^ + ^ = {Avx - A*v'^) C, which gives three more relations 





= AC, 


Ac) , Ab) 


= -A*C, 




= 0. 



(3.8) 



These relations express the Q s in terms of the Q s. Combining Eqs. ( p.7|) with 
yields 

Ac)_ wA*-A (,) 

d^^ = ^^^C-(--l)d^^ (3.9) 
w + 1 

At zero recoil, only d'^'' can give a non-vanishing contribution to the matrix elements of the 
weak currents in Eq. (|3.1|) . It is determined in terms of A* — A and ({1), since Eqs. ( p.8|) 
and (|3.9|) imply that 
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cr(i) = -cf(i) = ^c(i). 



(3.10) 



We use Eq. (|3.9|) to eliminate ^2'^'' and (^'^^ in favor of (i^^ and (. 

There are also order Aqcd/'^q corrections to the effective Lagrangian, given in Eq. ( [1.4D . 
Time ordered products involving Okin can be parameterized as 



d'x T {0Si,,(x) [h^:) r /^f ] (0) } = xS Tr {k^' T H^} . 



(3.11) 



These corrections do not contribute at zero recoil. The chromomagnetic corrections have 
the form 

d'x T {oLtg,,, (x) [h^} r h^:^] (0) } = Tr {n% K^. la^^ V 
^ / d'xT {Oltl^^ix) [h^^ r /.W] (0) } = Tr {u^ K^, T ^a-^H^], (3.12) 



In this case the most general form of TZ^^^ 

-oW (c) , (c) (b) , (6) / 



(3.13) 



At zero recoil the contribution of x^^^ vanish because fa(l + ^)cr"^(l + |^S) =0, while that of 
Xi'^^ vanish because (1- |(5)7„7^(1 + |^) = (1 - '^){-faVf3 - 7/3^'«)(l + ?^)- 

Using Eqs. ( |3.5| )-( prT2D , it is straightforward to express the form factors Qi parameterizing 
B —>■ D^eUe and B D^eVe semileptonic decays in terms of Isgur-Wise functions. The 
order Aqqd /"m-b Lagrangian corrections arise only in the combination Xb = Xkc + ^xf^ ~ 
2{w + l)x2^''- Dropping the c superscript from ([^'^ and Xi^\ "we obtain 



9+ = £c 



2(t/7-l)Ci-3C' 



wA* - A 
w + l 



A*(2w + 1) - A(w + 2) 



w + l 



C-2(w7-l)Ci 



9- = C + £c [xke + 6xi - 2(w + 1)^2] +ehXb- 
The analogous formulae for B ^ Dleue are 



(3.14) 



gA = ( + £c 



wA* -A 
w + l 



( + Xkc- 2xi 



Bb 



A*{2w + l)-A{w + 2) 
w + l 



C-2{w-l)Ci-Xb 



gv, = iw-l)C + e, [{wA* - A)C + iw- l)(xke - 2xi)_ 

-St {[A*{2W + 1) - A{W + 2)] C - 2{W^ _ 1) - («; - 1);^,} 

9V2 = 2ec (Ci - X2) , 

wA* - A 



9V3 = -C ~£ 
+£b 



w + l 

A*{2w + l)-A{w + 2) 



w + l 



C + 2Ci + Xke - 2xi + 2x2 

C,-2{w-l)C,-Xb 



(3.15) 
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These equations show that at zero recoil the leading contributions to gvi and (7+ of order 
A.QCJ) /tuq are determined in terms of A* — A and C(l)- Explicitly, 



-{e, + e,) (A*-A)C(1) 



9+{^) - 

gy^{l) = {e,-3e,) (A*-A)C(1). 



(3.16) 



For approximation A we shall again expand the double differential decay rates in Eq. ( p?2| 
in powers of w — 1, 



d^r 



D* 



dw dcos 9 



dw dcos 9 



3To r*3 (^2 _ ^)3/2 gi^2 ^ ^ (y; _ 1)- 4' 

n 

3To C'(l) rf V^I^^ E - 1)" i sin' 4"^ 



n) 



(3.17) 



+ (1 - 2r> + rf ) [(1 + cos^ 
The coefficients for the decay rate into Dq are 



4 cos 



SO) 



1 - ro)*' 1 + 2£,(xke + 6x1 - 4X2) + 45^6 + 3(£e + ^fc) (1 " r^') (A* - A) + . . . 



2(1 



„*\2 /■' 



C' + ... . 



(3.18) 



For the decay into Dl the coefficients are 



■'I 


= (ee - 3£b)2 (1 - 


r*)2 (A* - A)2 + . . . ^ 


,(1) 


= -2(£, - 3£fe) (1 


- rf) (A* - A) + ... , 


,(2) 


= (l + r*)2 + ... 






= (e, - 3€bY (A* 


-A)2 + ..., 




= 2 + 4(£e - SSb) 


(A* - A) + 4£e(Xke - 2Xl) + 4£bX6 + 


.(2) 


= 2(1 + 2C) + . . 






= (sc - 3eb) (A* - 


- A) + . . . , 




= 1 + ... . 





(3.19) 

Note that at zero recoil and at order w — 1 the contributions to decay proportional to 
A* — A depend on the anomalously small combination ec — <iSb ~ 0.05GeV~^. Thus Aqcd/'"^q 
corrections enhance B hj a much smaller amount than they enhance B ^ Di decay. 

On the other hand, the B Dq decay rate receives a large enhancement from Aqcb/it^q 
corrections, similar to B ^ Di. 

In approximation A, B ^ Dl is treated the same way as B ^ Di in Sec. II. B —>■ Dq 
is treated a.s B ^ D2 in Sec. II, since these rates contain an additional factor of — 1. 
Approximation B is also very similar to that in Sec. II, except that in the present case there is 
only one unknown Aqcd / f^Q Isgur-Wise function, (i (once time ordered products involving 
the chromomagnetic operator are neglected, and the matrix elements of the time ordered 
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products involving the kinetic energy operator are absorbed into the mg oo Isgur-Wise 
function, (). In approximation Bi we set Ci = in Eqs. ( p. 141) and ( 3.15 ). This is identical 
to saturating the first relation in Eq. by ([''\ i.e., setting ([^^ = A^. In approximation 



B2 we set Ci = in Eqs. ( p. 14] ) and ( |3.15| ), which is identical to setting ([''^ = 0. To the 
extent the first relation in Eq. (|3.8|) can be taken as a hint to the sign of (i, the difference 
between approximations Bi and B2 gives a crude estimate of the uncertainty related to the 
unknown Aqcb/i^q corrections. 

As in the previous section, the expression for the decay rate in terms of form factors 
in Eq. ( |3.2| ) implies that s^i^ /t^i^ = (1 — r^)^ to all orders in the Aqcd/''^q expansion. 
Thus the ratio of helicity zero and helicity one B Dl decay rates at zero recoil is 



lim„ 



(drgr°Vd^)/(dr^r'Vd^) = 1/2. 



(|A| = 



^1) 



A. Predictions 

A model independent prediction similar to that in Sec. II can be made for the slope 
parameter of semileptonic B decay into the helicity zero Dl. We write the semileptonic 
decay rate into the helicity zero Dl as 

— ^ = 4ro rf (1 - riy V^i;^C\l) (e, - 3e.)' (A* - A)^ [i _ p^, _ i) + . . . . 
aw L 1 J 

(3.20) 

The relationship between and s^i^ implies that the slope parameter pj^* for helicity zero 
Dl is 

1 -I- r* 9 

pI, = L±Il £^ ^ + 0(1) . (3.21) 

1 _ rl {e, - Ssb) (A* - A) ^ ^ ^ ^ 



As in Sec. II, this slope parameter is of order rriQ /Aqcd- It would be very hard experimentally 
to test this model independent prediction, since the Dl is expected to be of order 100 MeV 
broad, and also because Ec — 36^ is so small. 

Predictions for the B — > Dq e Ve and B ^ DleUg rates are shown in the first two columns 
of Table normalized to C^(l) times the measured B —>■ DieUe rate. These results are 
obtained using (' = —1, and A* — A ^ 0.35 GeV corresponding to 1 < w < 1.33. This 
value of A* — A has at least a 50 MeV uncertainty at present, as it follows from model 
predictions for the masses of the sj"^ = charmed mesons, fn^ ~ 2.40 GeV p^, and 



from the fact that = A'^ in nonrelativistic quark models with spin-orbit independent 
potentials. Although the Dl state is expected to be somewhat heavier than the Dq, we use 
the kinematic range 1 < w < 1.33 for both decays. The results in the first two columns of 
Table |V| are quite sensitive to the value of (' and Ci- In approximation Bi, for example, 
B{B ^ D*e z/e)/[C^(l) X 0.006] changes from 1.92 at C' = to 0.54 at C' = -2. In the same 
range of C', B{B D* e z/e)/[C^(l) x 0.006[ changes from 0.72 to 0.24. The effect of Ci is 
also important; in the range —0.75 GeV < < 0.75 GeV, the Dq and Dl branching ratios 
change from 1.68 to 0.66 and 0.30 to 0.63, respectively. Therefore, even if ( were known 
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Approximation 


X 0.006 


B{B Dleue) 
C^(l) X 0.006 






0.30 


0.66 


1.07 


Boo 


0.33 


0.46 


1.61 


A 


1.03 


0.65 


0.80 


Bi 


1.11 


0.44 


1.03 


B2 


0.85 


0.53 


1.05 



TABLE IV. The first two columns show semileptonic B branching ratios into and D\ 
normahzed to C^(l) times the measured branching ratio B{B Dievg) = 0.6%, assuming 
C' = C'(l)/C(l) = ~1- The sum of Dq + Dl rates relative to i? — > Di is in the third column, 
using the nonrelativistic constituent quark model prediction in Eq. (|3.22| ) and f' = —1.5. 



from models or lattice calculations, there would still be a factor of two uncertainty in the 
theoretical predictions for the semileptonic B —>■ Dq and rates; but the uncertainty in 
the sum of these two rates is smaller. 

To obtain even a crude absolute prediction for the B ^ Dl, Dq rates, a relation between 
the sy = and |^ Isgur-Wise functions is needed. In any nonrelativistic constituent quark 



model with spin-orbit independent potential, C, and r are related by pi, 19 



V3 



t{w) 



(3.22) 



since both of these spin symmetry doublets correspond to L 
imphes 



1 orbital excitations. This 



C = - + f'. 

^ 2 



(3.23) 



I doublet 



In the same approximation, f/^e = Xke-Q 

Predictions for the B semileptonic decay rate into the states in the 
that follow from Eq. (|3^) are shown in the last column of Table (For this quantity, 
approximations Bj (i = 1, 2) contain a somewhat ad hoc input of combining the Bj prediction 
in Sec. II with the Bj prediction for B Dq, Dl.) For f = —1.5, the ^ doublet contributes 
about 1.0 X B{B — > DieUe) ~ 0.6% to the total B decay rate. Varying ri2 and Ci in 
approximation B results in the range (0.6 — 1.7) x B{B DicVe) for the sum of the Dq 
and Dl rates. This combined with our results for R = T £,* /T in Sec. II is consistent with 



the ALEPH measurement of the branching ratio for the sum of all semileptonic decays 
containing a D^*'' vr in the final state to be (2.26 ± 0.44)%. 

The semileptonic decay rate into D and D* is about 6.6% of the total B decay rate 0. 
Our results then suggest that the six lightest charmed mesons contribute about 8.2% of the 
B decay rate. Therefore, semileptonic decays into higher excited states and non-resonant 



^A relation between ri^2 ^-ncl C,i may also hold in this model. 
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multi-body channels should be at least 2% of the B decay rate, and probably around 3% 
if the semileptonic B branching ratio is closer to the LEP result of about 11.5%. Such a 
sizable contribution to the semileptonic rate from higher mass excited charmed mesons and 
non-resonant modes would soften the lepton spectrum, and may make the agreement with 
data on the inclusive lepton spectrum worse. Of course, the decay rates to the broad I"*" 
states would change substantially if the nonrelativistic quark model prediction in Eq. (|3.22|) 
is wrong. Semileptonic B decay rate to the six lightest charmed mesons could add up to 
close to 10% if ( were enhanced by a factor of two compared to the prediction of Eq. ( p. 221) . 
However, model calculations seem to obtain a suppression rather than an enhancement 
of ( compared to Eq. ( |3.22| ). Thus, taking the measurements for the B —>■ D, D*, and 
Di semileptonic branching ratios on face value, a decomposition of the semileptonic rate as 
a sum of exclusive channels seems problematic both in light of our results and the above 
ALEPH measurement for the sum of all semileptonic decays containing a D^*^ vr in the final 
state. 

IV. OTHER EXCITED CHARMED MESONS AT ZERO RECOIL 

In the previous two sections matrix elements of the weak vector current and axial- vector 
current between a B meson and an excited charmed mesons with sj' = l"*^ and quantum 
numbers were considered. Here we consider such matrix elements at zero recoil for excited 
charmed mesons with other s^^ quantum numbers. Only charmed mesons with spin zero or 
spin one can contribute at this kinematic point. The polarization tensor of a spin n state 
is rank n, traceless and symmetric in its indices, and vanishes if it is contracted with the 4- 
velocity of the state. For matrix elements of the axial- vector or vector current, at least n — 1 
indices of the charmed meson polarization tensor are contracted with v'^, the four velocity 
of the B meson. Consequently, for n > 1 these matrix elements vanish at zero recoil, where 
V = v' . In this section we work in the rest frame, v = v' = (1,0), and four-velocity labels 
on the fields and states are suppressed. 

For spin zero and spin one excited charmed mesons, the possible spin parities for the 
light degrees of freedom are = which we have already considered in the previous 

sections, and s^^ = | , | . In the nonrelativistic constituent quark model, the | states 
are interpreted as radial excitations of the ground state {D, D*) doublet and the | states 
are L = 2 orbital excitations. In the quark model, these states are typically expected to 
be broad. The mass of the lightest sf' = | doublet is expected around 2.8 GeV, while the 
lightest excited states with sj'- = \ are expected around 2.6 GeV [^.0 {B decays into 
radial excitations of the s^' 7^ | states have similar properties as the decay into the lightest 
state with the same quantum numbers.) 



®The lightest | states may be narrow since decays to the sj^ = ^ and | multiplets are 
suppressed by the available phase space, and decays to D*^*-' tt in an 5-wave are forbidden by 
parity. 
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In the TTiQ oo limit, the zero recoil matrix elements vanish by heavy quark symmetry. 
For the excited s^* = | states, the mg oo Isgur-Wise functions vanish at zero recoil due 

to the orthogonality of the states. The matrix elements for the sj* 7^ | states vanish at zero 
recoil due to spin symmetry alone, and therefore the corresponding mg — > 00 Isgur-Wise 
functions need not vanish at zero recoil. 

Using the same methods as in Sections II and III, it is straightforward to show that 
Aqcd/'^q corrections to the current do not contribute at zero recoil. For the s]' = | 
states, this follows from the heavy quark equation of motion. For the s'^' = | states, the 
Aqcb/^q corrections to the current can be parameterized similar to Eqs. (|2.^ ) and ( p.9|) . In 
this case the analogue of in Eq. ( p.5| ) satisfies = = ^. Recall that the t^'^^Qctx 
in Eq. (|2.9| ) was the only term whose contribution at zero recoil did not vanish due to the 
Vf^F^ = property of the Rarita-Schwinger spinors. Here, the analogous term is placed 
between 1 — and 1 + and therefore also disappears at f = v'. 

It remains to consider the Aqcb/^q contributions to the | and | matrix elements 
coming from corrections to the Lagrangian in Eq. ( p.. 41) . These are written as time ordered 
products of Ol^^{x) and O^Jg{x) with the leading order mq — > 00 currents (e.g., Eq (|2.16|) ). 
At zero recoil it is useful to insert a complete set of states between these operators. Since 
the zero recoil weak currents are charge densities of heavy quark spin-flavor symmetry, only 
one state from this sum contributes. For the sj*^ = | multiplet this procedure gives 

{D<-){e)\A\B) -e fM + 3 \ {D<-Ke)\ 0^1,(0) \D*ie)) 



y/rrij^^Cn) niB (A(") — A) [ \2mc 2m;,/ rriD 



1 \ {D*'^-\^)\Ol^liO)\D*i^)) \ .4^^. 

\2mc 2mi)J y/mjj,(n) mc* J 



and 



{D^-)\V^\B) 1 / 1 , 1 \ (o) + o£i(0)|D) 



iym£)(n) itlb (A(") — A) V 2mc 2mbJ ^m£){„) mr, 



(4.2) 



Here we have denoted spin zero and spin one members of the excited sY = \ multiplet 
by D^"-' and respectively, and the analogues of A by A^"^ Heavy quark spin-flavor 

symmetry was used to write the effects of O^;^ and O^^g in terms of matrix elements of O^j^ 
and O^^g. This neglects the weak logarithmic dependence on the heavy quark mass in the 
matrix elements of Omag- For the spin one member of the = | multiplet, which we 
denote by 



{Drie)\A\B) _ -e f 1 ^ {Dl*{s)\ 0(^(0) \D*{e)l ^^^^ 



^mDi* ms {A** - A) \2mcJ ^yrnj^f^lnj^ 

For the sj' = | and | excited charmed mesons, the correction to the Lagrangian, SC 
in Eq. ( |1.4| ), gives rise to an order Aqcd/'^tIc contribution to the matrix elements of the weak 
currents at zero recoil. Formulae similar to those in Eqs. ( |4.1| )-( |0| ) hold in the sj* = 
1^ cases, but the corresponding matrix elements vanish due to the parity invariance of the 
strong interaction. 
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Approximation 




B{B Ditt) 


B{B 
B{B 


^ D2 vr) 






0.39 




0.36 






0.26 




1.00 


A 




0.29 




0.21 


Bi 




0.19 




0.41 


B2 




0.20 




0.56 


TABLE V. 


Predictions 


for the ratios of branching 


ratios, B{B Di-k)/B{B DieUe) and 



B{B D2'7t)/B{B Diir), using factorization and assuming f' = r'(l)/r(l) = —1.5. 



V. APPLICATIONS 
A. Factorization 

Factorization should be a good approximation for B decay into charmed mesons and a 
charged pion. Contributions that violate factorization are suppressed by Aqcd divided by 
the energy of the pion in the B rest frame or by asimq). Furthermore for these decays, 



factorization also holds in the limit of large number of colors. Neglecting the pion mass, the 
two-body decay rate, F^, is related to the differential decay rate dTgi/dw at maximal recoil 
for the analogous semileptonic decay (with the vr replaced by the eue pair). This relation is 
independent of the identity of the charmed meson in the final state, 

mfj r \ aw J 

Here r is the mass of the charmed meson divided by itlb, u^max = (1 + '"^)/(2r), and /^r ~ 
132 MeV is the pion decay constant. C is a combination of Wilson coefficients of four-quark 
operators, and numerically C \ Vud\ is very close to unity. 

These nonleptonic decay rates can therefore be predicted from a measurement of dFsi/dty 
at maximal recoil. The semileptonic decay rate near maximal recoil is only measured for 
B D^*^ eUe at present. The measured B D^*'> tx rate is consistent with Eq. ( p.l[ ) at 
the level of the 10% experimental uncertainties. In the absence of a measurement of the 
B — > {Di.DD eVf. differential decay rates, we can use our results for the shape of dFgi/dtu 
to predict the B ^ D\ti and B D\t[ decay rates. These predictions depend on the 
semileptonic differential decay rates at Wmax, where we are the least confident that Aqcd/'^q 
terms involving A and A' are the most important. With this caveat in mind, we find the 
results shown in Table 0. 

At present there are only crude measurements of the BiB —>■ Di n) and B{B tt) 
branching ratios. Assuming -B(Di(2420)° ^ D*+ tt") = 2/3 and B{D*{2AQ0f D*+ tt") = 
0.2, the measured rates are 

B{B- Di(2420)°7r-) = (1.17 ±0.29) x 10-^ 

B{B- D;{2A60f 71-) = (2.1 ± 0.9) x 10"^ . (5.2) 
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FIG. 5. Factorization prediction for B{B Ditt)/B{B Dieve) as a function of 
l)/r(l). The dotted curve shows the mq oo hmit (Bqo), sohd curve is approxima- 
tion Bi, dashed curve is B2. 

A reduction of the experimental uncertainty in B{B D2 vr) is needed to test the prediction 
in the second column of Table |V|. 

The prediction for B{B Di7i)/B{B DicVf.) in approximation B is fairly in- 
dependent of Ti 2, but more sensitive to f'. The latter dependence is plotted in Fig. |^ 
for > f' > —2. Not absorbing ?7ke into r results in the following weak dependence: 
B{B Di n)/B{B ^ Z^i e z/e) oc 1 + 0.27 ec%c + • • •• Assuming that factorization works at 
the 10% level, a precise measurement of the B{B Di it) rate may provide a determination 
of f'. The present experimental data, B{B Diti)/B{B DieVe) — 0.2, does in fact 
support f' ~ —1.5, which we took as the "central value" in this paper, motivated by model 
calculations. 

The prediction for B{B D2 'k)/B{B -Di tt), on the other hand, only weakly depends 
on f', but it is more sensitive to ti^2- Varying ri^2 in the range —0.75 GeV < fi_2 < 0.75 GeV, 
we can accommodate almost any value of B{B D2'n)/B{B Din) between and 
1.5. This quantity depends more sensitively on ti than on T2. In Fig. |^ we plot B{B — 
D2 7i)/B{B DiTi) in approximation B as a function of fi setting 7^2 = (solid curve), 
and as a function of f2 setting fi = (dashed curve). Not absorbing ?7ke into r results in 
the following dependence: B{B DIti)/B{B — > Di tt) oc 1 + 0.75ec%e + • • •• This ratio 
and R depend on r)ke and fi. In the future experimental data on these ratios may lead to a 
determination of i7ke and fi. 

If the experimental central value on B{B —>■ D2 tt) does not decrease compared to 
Eq. ( |5.2| ), then it would suggest a huge value for fi, leading to a violation of the ALEPH 
bound on R (see Fig. 0). The approximation B results in Tables ^ and |V| can be com- 
bined to give B{B D2n)/B{B D2eUe) = 0.15. Varying fj, r/ke and f does not bring 
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FIG. 6. Factorization prediction for B{B — > D2Tt)/13{B Ditt) as a function of fi(= ri/r) 
for f2 = (solid curve), and as a function of f2 for fi = (dashed curve). 

this quantity close to the current experimental limit. Therefore, if the branching ratio for 
B e is below the ALEPH bound, then B{B D2 vr) should be smaller than the 

central value in Eq. (|5.2| ). 



B. Sum Rules 

Our results are important for sum rules that relate inclusive B ^ XcCUe decays to the 
sum of exclusive channels. The Bjorken sum rule bounds the slope of the B D^*^ e Ue 
Isgur-Wise function, defined by the expansion ^{w) = 1 — {w — 1) + . . . . Knowing 
would reduce the uncertainty in the determination of \Vcb\ from the extrapolation of the 
B D^*'> e Vf. spectrum to zero recoil. The Bjorken sum rule p9|,p!9 



IS 



p^ = i + E^4^ + 2E^^^ + ---- (5.3) 

Throughout this section the ellipses denote contributions from non-resonant channels. C^^^ 
and t'^p^ are the Isgur-Wise functions for the exited sY = I"*" and |^ states, respectively (for 
m = p = these are the orbitally excited states discussed in Sec. II and III, and m,p > 1 
are radial excitations of these). Since all terms in the sums, as well as the contributions 
replaced by eUipses, are non-negative, a lower bound on can be obtained by keeping only 



Hn Ref. ig |C(""ni)lV4 was denoted by ItI"'^ , and |r(p)(l)|V3 was denoted by |r3^^J(l)p. 
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the first few terms on tlie riglit-fiand-side of Eq. ( |5.3| ). Using Eqs. ( p.31| ) and ( p.22|) , we find 

and I 



tliat tlie contribution of tlie lowest lying sj' = and |^ states implies the bound 



,.>i + Ml)£ + 2Mi)£^0,75. (5.4) 



4 4 



The contribution of the states through ({1) to this bound, which rehes on the quark 
model result in Eq. ( p.22|) , is only 0.17. 



An upper bound on p follows from an upper bound on the excited states contribution 



to the right-hand-side of Eq. (p.3|). This sum rule was first derived by Voloshin |30 



-A = 5](A*(™) -A)^ — ^ + 2^(A^(P) - A) ' +... . (5.5) 

Here A*^'") and A'^^^ are the analogues of A* and A' for the exited s'^' = |^ and |^ states, 
respectively. Eq. ( |5.5| ) combined with Eq. ( pl3| ) implies that < 1/4 + A/(2£:i), where ei is 
the excitation energy of the lightest excited charmed meson state. However, knowing ^(1) 
and r(l) does not strengthen this bound on significantly. On the other hand, Eq. ( p.5[ ) 
imphes the bound A > 0.38 GeV (neglecting perturbative QCD corrections). The model 
dependent contribution of the states to this bound is only 0.12 GeV; while the bound 
A > 0.26 GeV from only the |^ states is fairly model independent. 

A class of zero recoil sum rules were considered in Ref. |^T|. The axial sum rule, which 



bounds the B — > D* form factor (that is used to determine |Vcf,|) only receives contributions 
from = I and | states, which were discussed in Sec. IV. It has the form 



|2 



\{X,ie)\A\B)\^ _ 2 As A1 + 3A2 / 1 1 



^B^D'{^)\ +J2 J =Va 2+ ~A \~ + ~ + ^ ' (^-6) 



where rjA is the perturbative matching coefficient of the full QCD axial-vector current onto 
the HQET current, denotes spin one states (continuum or resonant) with s^*^ = ^ and 

I , and Fb_>z)*(1) is defined by 

{D*{e)\A\B) ^ 



[l)e. (5.7) 



Neglecting the contributions of the excited states Xc to the left-hand-side, gives an upper 
bound on \Fb^d*{^)\'^- Using the nonrelativistic constituent quark model, we estimate using 



Eq. (O) that the contribution of the first radial excitation of the D* to the sum over X^ in 
Eq. is about 0.1. For this estimate we took A^^) - A = 450 MeV, O^^l^ = C 6^{r) Sc ■ Sq 

(fixing the constant C by the measured D* — D mass splitting), O^;^ = V^, and used 



the harmonic oscillator quark model wave functions of Ref. p2|. A 0.1 correction would 



significantly strengthen the upper bound on Fb-*d*{^) and have important consequences for 
the extraction of the magnitude of Vcb from exclusive B D*eUe decay. Note that s^* = | 
states do not contribute to the zero recoil axial sum rule in the quark model, because their 
spatial wave functions vanish at the origin. 
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The = 1+ members of the s^^ = I"*" and sj*^ = I"*" doublets contribute to the vector 



sum rule, which is used to bound Ai. This sum rule reads 

\2 \/-(m)/T\\2 



4 f^c p 



4ml 



+ 



A2 Ai + 3A2 
2 ~ ' 



m 



11 2 
4 Vm^ Srricmb 



Ai = -3E(A 



* (m) 



A)^ 



|^M(l)p 



6^(A'{p) _ A) 



r(p)(l)p 



+ 



(5.8) 



This relation can be simplified by setting mb/rric to different values. Taking mt = rric yields 



(5.9) 



whereas rric^ rrib^ Aqcd gives 



Ai + 3A2 = -95] (A* ('"^ - A)' 



+ . . . . 



(5.10) 



These relations can be combined to obtain a sum rule for A2, 
A2 = -2 5] (A* ('"^ - A)2 ^ '^ + 2 ^ (A'(f) - A)2 



-{p)(l)p 



+ 



(5.11) 



Eqs. ( ^.9|) and (|5.11|) were previously obtained in Ref. ||3^ using different methods. The 
strongest constraint on Ai is given by Eq. ( ^.10| ) (the sum rule in Eq. (|5.9|) only implies 
— Ai > (0.06 + 0.15) GeV^). Including the contribution of the lightest = \^ doublet to 
Eq. ( FTOD yields 



Ai < -3A2 -9(A* - A) 



2lC(l)P 



~ -3A2 -0.18GeV^ 



(5.12) 



neglecting perturbative QCD corrections. Note that only the broad D\ state (and its radial 
excitations) contribute to this sum rule, so the result in Eq ( 5.12 ) is sensitive to the relation 
between r(l) and C(l) in Eq. ( ^1221) . 

Perturbative corrections to the sum rules in this section can be found in Ref. ||33|1. 



VI. SUMMARY AND CONCLUSIONS 

The branching ratios for B ^ D eVe and B D* e Ve are (1.8 ± 0.4)% and (4.6 ±0.3)%, 
respectively ||^. This implies that about 40% of semileptonic B decays are to ex- 
cited charmed mesons and non-resonant final states. An excited charmed meson doublet 
(Di (2420), 1^2(2460)) with sY = |^ has been observed. These states are narrow and have 
widths around 20 MeV. With some assumptions, the CLEO and ALEPH collaborations have 
measured about a (0.6 ±0.1)% branching ratio for B ^ Die Ve- The decay B D\ e has 
not been observed, and CLEO and ALEPH respectively report limits of 1% and 0.2% on 
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its branching ratio. A detailed experimental study of semileptonic B decays to these states 
should be possible in the future. 

The semileptonic B decay rate to an excited charmed meson is determined by the corre- 
sponding matrix elements of the weak axial- vector and vector currents. At zero recoil (where 
the final excited charmed meson is at rest in the rest frame of the initial B meson), these 
currents correspond to charges of the heavy quark spin-flavor symmetry. Consequently, in 
the TTiQ — > oo limit, the zero recoil matrix elements of the weak currents between a B meson 
and any excited charmed meson vanish. However, at order Aqcd/''^(3 these matrix elements 
are not necessarily zero. Since for B semileptonic decay to excited charmed mesons most 
of the available phase space is near zero recoil, the Aqcd/"t.q corrections can play a very 
important role. In this paper we studied the predictions of HQET for the B DieUe 
and B — > D2ei?e differential decay rates including the effects of Aqcd / n^-Q corrections to 
the matrix elements of the weak currents. Since the matrix elements of the weak currents 
between a B meson and any excited charmed meson can only be nonzero for spin zero or 
spin one charmed mesons at zero recoil, the Aqcb/ttiq corrections are more important for 
the spin one member of the sj^ = |^ doublet. 

The Aqcb/itt-q corrections to the matrix elements of the weak axial- vector and vector 
currents can be divided into two classes; corrections to the currents themselves and correc- 
tions to the states. For B semileptonic decays to the Di, parity invariance of the strong 
interactions forces the corrections to the states to vanish at zero recoil. Furthermore, the 
corrections to the current give a contribution which at zero recoil is expressible in terms of 
the leading, mg oo, Isgur-Wise function and known meson mass splittings. This cor- 
rection leads to an enhancement of the B semileptonic decay rate to the Di over that to 
the With some model dependent assumptions, we made predictions for the differential 
decay rates for B — > Di e Pg and B D2eUe and determined the zero recoil value of the 
leading mq — > oo Isgur-Wise function from the measured B to Di semileptonic decay rate. 
The influence of perturbative QCD corrections on these decay rates were also considered 
but these are quite small. 

Factorization was used to predict the rates for the nonleptonic decays B ^ Din and 
S — > D2 TT. The ALEPH limit on the semileptonic decay rate to D2 implies a small branching 
ratio for B ^ D^n. The ratio B{B — > Di n)/B{B Die Pe) can be used to determine f'. 
The present experimental value for this quantity favors f' near —1.5. 

The most significant uncertainty at order Aqcd /i^Q arises from fi and r)kc- It may be 
possible to determine these quantities from measurements oi R = rD*/^Di and B{B — > 
D2 n)/B{B -Di vr). The u^-dependence of the semileptonic decay rates can provide impor- 
tant similar information. 

A broad multiplet of excited charmed mesons with masses near those of the Di and D2 
is expected. It has spin of the light degrees of freedom = \ , giving spin zero and spin 
one states that are usually denoted by Dq and D^. We studied the predictions of HQET for 
the B ^ DleUf. and B ^ Dlev^ differential decay rates including the effects of Aqcd /mq 
corrections to the matrix elements of the weak current. The situation here is similar to 
that in the case of the sj^ = I"*" doublet. Using a relation between the leading, mq 00, 
Isgur-Wise functions for these two excited charmed meson doublets that is valid in the 
nonrelativistic constituent quark model with any spin-orbit independent potential (and a 
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Approximation R = Tj:,*iTj:,-^ t(1) 

B^^ L65 1.24 4.26 

Bi 0.52 0.71 2.55 

B2 0.67 0.75 2.71 



TABLE VI. Predictions for r^ij/r^,^, t(1), and Td^+d*+di+di/^Di using f = -1.5. The 
results in the last column assume the nonrelativistic quark model prediction in Eq. ( p.22| ). 



6.0 X 10"^ 
BiB Dieue) 



1/^ 



^Di+D*+Dl+D*/'^Di 



few other assumptions), we determined the rates for B semileptonic decays to these excited 
charmed mesons. We find that branching ratio for B semileptonic decays into the four states 
in the s^'^ = and |^ doublets is about 1.6%. Combining this with the measured rates 
to the ground state D and D* implies that more than 2% of the B meson decays must be 
semileptonic decays to higher mass excited charmed states or nonresonant modes. Some of 
the more important results in Tables |I| and are summarized in Table [V^. 

We considered the zero recoil matrix elements of the weak currents between a B meson 
and other excited charmed mesons at order Aqcd/'"^q- Only the corrections to the states 
contribute and these were expressed in terms of matrix elements of local operators. 

Our results have implications for B decay sum rules, where including the contributions of 
the excited charmed meson states strengthens the bounds on (the slope of the Isgur-Wise 
function for B — > D^*^ e z/g), on Ai, and on the zero recoil matrix element of the axial- vector 
current between B and D* mesons. The latter bound has implications for the extraction of 
\Vcb\ from exclusive B D* eue decay. 
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APPENDIX A: PERTURBATIVE ORDER a, CORRECTIONS 



In this Appendix we compute order as and order asAQCD/''^Q corrections to the B —>■ 
{Di, D2) e Ve form factors. At this order both the current in Eq. (|2.7| ) and the order Kop-QjmQ 
corrections to the Lagrangian in Eq. ( |1.4| ) receive corrections. Matrix elements of the kinetic 
energy operator, , enter proportional to r to all orders in due to reparameterization 
invariance 
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The matrix elements involving the chromomagnetic operator are probably 
very small and have been neglected. Order as corrections to the b ^ c flavor changing 
current in the effective theory introduce a set of new operators at each order in Aqcd/i^q, 
with the appropriate dimensions and quantum numbers. The Wilson coefficients for these 



operators are known w-dependent functions |T7| , |35[] , which we take from p6[ . 



The vector and axial- vector currents can be written at order a. as 
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2 TTLr 



+ 



2mh 



TT 



+ 



7 75 ^ K 



2mr 



2mb 



TT 



+ . . . 



(Al) 



where the elhpses denote terms higher order in and Aqcd/^^q- Superscripts (1) denote 
corrections proportional to a^, 



The terms with superscript (2) in Eq. 



^m(2) 



2mb 

2mc 

iA 
2 nib 

i% 

2mr 



75 ft;,"' • (A2) 
denote corrections proportional to as Aqcd/'o^q, 



7^ + 2v'^ 



/X A 



2cv^g^'^+[-^'^ + 2v 



dwj 

CviY + CV2V^ + CV3V''') 



(A3) 



dwJ 



(c^i7^ + CA2V^ + CAsv'") 75 ^7^ + 2t; 
2cA3g^^ + (7^ + 2t;^ ^) (cai7^ + 0^2^^^ + 



+ 2CA2(?^S5 
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In these expressions the covariant derivatives, Dx, act on the fields h^^^ or h^^) , and partial 
derivatives with respect to w, d/dw, act on the coefficient functions cyj(w) and CAi{w). 
Using Eqs. ( |A^ ) and ( [A3|) it is straightforward to include the order as and ag AqcB/mQ 
corrections using trace formalism presented in Sec. II. The corrections with superscript (1) 
simply change the form of F in Eq. (|2.6| ), while those with superscript (2) change F in 
Eq. (^. 

The B DieUe form factors were defined in Eq. (|2.1| ), and their expansions in terms 
of Isgur-Wise functions at leading order in as were given in Eq. ( |2.2CI| ). The order as and 
order Aqcd/"^q corrections modify the results for /j in Eq. (|2.2CI|) to /j + {as/n) 6 ft. The 
functions 6fi are given by 

VeSfA = -{w + l)cAir - 2e, {wK' - K)[2cai + {w + l)c'^i + caz]t 

+ ec{w - 1) |[3cai - 2{w - 1)ca3]ti - (3cai + 4cA3)r2} 

- Sb [(A' + A)(w - l)cAi - 2(A' - wK){w + + 2 {wK' - K)ca2]t 

+ eb{w- 1) {[{2w + l)cAi - 2(w; - l)cA2]ri + (c^i - 40^2)^2}, (A4) 

Ve^/y^ = (1 - w^)cyir - 2ec {wA' - K){w + 1) [2cyi + {w- l)c!y^ + 2cy3]r 

+ Ec {W^ - l){[3cyi + 2{w + 2)cy3]ri - (3cyi + 2cy3)r2} 

-eb{w + 1) [(A' + K){w-l)cvi- 2(A' - wK) {w - l)c'yi + A{wA! - K)cv2] r 

+ Eb {W^ - 1){[{2W + l)cyi + 2(W + 2)cy2]ri + (Cyi - 2cy2)r2}, (A5) 
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V6Sfv, = -[3cvi + 2{w + l)cv2]r - 2e^ {wK' - A) \^dy^ + 2cv2 + 2{w + l)c'^2j 

- e,[[{Aw - l)cvi - 2{2w + l){w - l)cv2 - 2{w + 2)cv3]ri 

+ [5cvi + 2{l-w)cv2 + 2cv3]r2} 

- Ek {3(A' + A)cyi - 6(A' - ^A)cVi + 2[{w - 1)A' + (3w + l)A]cy2 

- A{K' - wK){w + l)c'y^Y 
+ eb[[?,{2w + l)cyi + 2{2w^ + l)cv2]ri + [3cyi + 2{w - 2)cv2]t2], (A6) 

VQ5fv, = [{w - 2)cyi - 2{w + l)cv3]r 

+ 2£c (wA' - A){2cyi + {w- 2)cyi - 2[cv3 + iw + 1)^3] }r 

+ £^{[(2 + w)cvi + 2{w^ -3w- l)cv3]ri + [{Sw + 2)cyi + {Aw - 2)cv3]r2} 

+ €b [(A' + A){w + 2)cvi + 2(A' - wA)i2 - w)c'y^ + AA'{w + l)cv2 

- 2{A' + A){w- l)cv3 + 4(A' - wA){w + 1)0^3] r 
-6b{ [i2w^ + 5w + 2)cvi + 2«;(2 + w)cv2 + 2{l+w- 2w'^)cv3\ri 

+ [(2 + w)cvi - 2wcv2 - 2{w - I)cv3]r2}- (A7) 

Here cvi and cai are functions of w, and prime denotes a derivative with respect to w. Note 
that at zero recoil is known in terms of A' — A and r(l), as expected from our results 
in Sec. II. 

For B D^eVe decay, the and order ag AQCj^/mq corrections modify the leading 
order form factors in Eq. ( p.21|) to ki ki + (as/vr) 5ki. The functions 5ki are 



5k\ 



-CviT - 



2cy^{wA' - A)r + (cvi - 2wcv3)ri - {cyi + 2cv3)r2 



(A8) 



£b 



{[(A' + A)cvi - 2(A' - wA)c'y^]T - [{2w + l)cvi + 2wcv2]ri - (cyi + 20^2)^2}, 



Sk, 



(W + 1)caiT - Be 



2(c'^i + wd^i - ca3){wA' - A)r + {w - l)cAi(ri - T2] 
eb{[{A' ^A){w- l)cAi - 2(A' - wA){w + l)c'^i 



+2{w' - l)cA3ri 

-2(wA' - A)cA2]r - - l)[cAi(ri + ts) + 2{wcai - wca2 - CA2)ri] k (A9) 



5kA2 = CA2T + Be {2c^2(^A' - A)t - [2cai - {2w + 1)ca2 + 2cA3]ri + CA27-2} 

+ eb{[{A' + 3A)cA2 - 2(A' - wA)c'^^]t - {2w + 3)0^2^1 - CA2T2}, (AlO) 

6kA3 = {CAI + CA3)t + Sc 2(c'^i + 0^3) (w A' - A)r - (c^l - CA3)(t"1 + T2) + 

+ eb{[{A' + A){cAi + CA3) - 2A'cA2 - 2(A' - wA){c'^, + c'a,)]t 

- {cai + CA3){ri + T2) - 2w{cai - Ca2 + cas)^}- (All) 
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To compute the corrections to the results obtained in Sec. II, it is sufficient to expand 
the Wilson coefficients cyi and cai to linear order in w. We take cyi and caz and their ffist 
derivatives at zero recoil from Ref. |3^. To evaluate these, we choose to integrate out the c 
and b quarks at a common scale /i = y/m^Pmb, giving for Cyi and Cai 

cyi(l) = -^-^ln^~0.91, 
2(1 - z + z\nz) 



2z(l - z + lnz) 



CAi(l) = -^-[^lnz~-0.42, 

_ 2[3~2z-z^ + {5-z)zlnz] ^ 

CA2{1)- 3.^_ .3 - 1.20. 



2z\l + 2z-3z^ + (5z-l)\nz] , 
CAsil) = ^ - 0.42 . (A12) 



The derivatives Cy^ and c^j at zero recoil are 

, 2[13-9z + 9z^ -13z'^ + 3{2 + 3z + 3z^ + 2z^)\nz] ^ 
cv.il) = 0.20 , 

2[2 + 3z-6z' + z' + 6z Inz] ^ 



, 2z[l-6z + 3z^ + 2z^ -6zHnz] ^ 

C,/q(l) = ; ~ 0.05 , 

^3V ; 9{l-z)^ 

, 2[7 + 9z-9z^-7z^ + 3{2 + 3z + 3z^ + 2z'')\nz] 
CMil) = ^YJT^^ " 



, 2\2-33z + 9z'^ + 25z^ -3z'^ -6z(l + 7z) \nz] ^ ^„ 
CA.il) = ^j^, 0.37, 

c',3(l) - _2-[3-25.-9.- + 33.3-y-6.V + .) m.] ^ _^ ^^^^^ 

9(1 z) 

Here z = rric/mi,, and the numbers quoted are for z = 1.4/4.8. 

Using these values and the as corrections for the form factors above, we find the correc- 
tions given in Table ^ to the leading order results summarized in Table fT[ 
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